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WE: 离子 束 放 射 治疗 (放疗 ) 越 来 越 广泛 地 被 运用 于 肿瘤 治疗 ， 目 前 使 用 传统 X 射线 计 
算 机 断层 扫描 《CT) 来 制定 离子 放疗 治疗 计划 时 存在 相对 阻止 本 领 (RSP) 具有 较 大 误差 
的 问题 ， 理 想 的 方式 是 直接 使 用 高 能 离子 束 对 患者 成 像 ， 生 成 离子 放疗 治疗 计划 设计 时 所 
需 的 医学 图 像 ， 避 免 RSP 转换 从 而 减 小 RSP 的 误差 。 本 研究 使 用 蒙特 卡 罗 程 序 
Geant4/Gate 搭建 了 离子 CT 模拟 平台 ,设计 了 理想 与 现实 两 种 离子 CT 成 像 系 统 探测 器 ， 
利用 最 大 似 然 法 和 ASD-POCS 算法 实现 了 图 像 重建 ， 研 究 了 理想 与 现实 两 种 设置 、 多 种 和 
量 以 及 不 同 离子 种 类 对 模 体内 硫 、 磷 、 碳 和 钙 四 种 材料 插件 重建 RSP 误差 的 影响 。 结 果 表 
A, 330MeV 质子 在 理想 与 现实 设置 两 种 情况 下 的 RSP 相对 误差 都 小 于 1.547%， 理 想 设置 
下 的 RSP 相对 误差 远 小 于 现实 设置 ， 现 实 设置 下 各 材料 的 RSP 重建 误差 接近 理想 设置 下 
的 三 倍 ; 质子 的 RSP 相对 误差 随 着 能 量 的 增加 而 减 小 ， 在 230MeV 时 的 RSP 相对 误差 最 
大 ， 在 硫 、 磷 、 碳 和 和 钉 四 种 材料 中 分 别 为 2.855%、2.468%、1.653%、2.553%。 在 330MeV 
时 的 RSP 相对 误差 达到 最 小 ， 在 硫 、 磷 、 碳 和 钙 四 种 材料 中 分 别 为 0.181%、0.351%、 
0.250% 和 0.245%; 在 能 量 330MeV/u 下 ， 碳 离子 在 硫 、 磷 、 碳 和 钙 四 种 材料 中 的 RSP 相对 
误差 分 别 为 0.060%、0.281%、0.150%、0.082%， 误 差 均 在 0.281% 以 内 ，RSP 相对 误差 小 
于 330MeV 的 质子 。 因 此 ， 与 质子 CT 相 比 ， 碳 离子 CT 能 为 离子 放疗 治疗 计划 设计 提供 
更 为 精确 的 RSP 数据 。 
关键 词 ， 离子 CT; 蒙特 卡 罗 模 拟 ， 质 子 ; 碳 离子 ， 相 对 阻止 本 领 
中 图 分 类 号 : R814.42 文献 标志 码 : A doi: 


1 引言 


离子 束 〈 如 质子 、 碳 离子 ) 放射 治疗 是 治疗 深层 肿瘤 的 一 种 有 效 方法 。 相 比 传统 用 光 
子 来 进行 放射 治疗 ,离子 束 〈 如 质子 、 碳 离子 ) 治疗 改善 了 对 肿瘤 的 治疗 ,健康 组 织 受 到 的 辐 
射 更 少 ,造成 的 放疗 副作用 也 越 小 。 这 是 因为 与 传统 的 光子 放射 治疗 ， 离 子 束 放射 治疗 拥有 
独特 的 物理 和 放射 生物 学 原理 来 帮助 提高 治疗 效果 。 离 子 束 治疗 是 一 种 新 型 放射 治疗 手 
段 ， 具 有 独特 的 物理 剂量 -深度 曲线 ， 即 布拉格 峰 ， 它 的 横向 散射 和 射程 歧 离 小 ， 在 其 射程 
末端 沉积 大 部 分 能 量 ， 具 有 人 尖锐 的 高 剂量 布拉格 峰 ， 可 以 精确 进行 肿瘤 靶 向 治疗 ， 从 而 拥 
有 更 高 的 均匀 剂量 ， 同 时 极 大 减少 周围 正常 组 织 的 损伤 ， 呈 现 较 高 的 生物 学 效应 。 据 
PTCOG (国际 离子 治疗 联合 会 ) 发 布 的 最 新 数据 ， 截 至 2021 年 底 ， 全 球 共 有 32 万 多 例 患 
接受 了 离子 治疗 ， 其 中 接受 质子 治疗 的 患者 28 万 多 例 ， 占 总 数 的 86.1%; 接受 碳 离子 治疗 
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的 患者 4 万 多 例 ， 占 总 数 的 12.9%; 其 余 3500 患者 接受 了 和 氨 离 子 、 
中 。 离 子 治疗 相 比 传统 的 光子 治疗 主要 有 三 个 优势 QB1。 
的 剂量 较 低 ， 不 仅 能 有 效 减少 短期 的 不 良 
者 生活 质量 ， 减 少 死亡 率 ， 5 
限度 的 投放 到 肿瘤 上 ， 改 善 肿 瘤 治 疗 。 
的 时 间 内 更 有 效 的 消灭 目标 内 的 肿 ; 


副 作 月 


介子 和 其 他 离子 治疗 


第 一 个 优势 是 离子 治疗 对 健康 组 织 
昌 ， 在 完成 治疗 后 数 年 内 ， 也 能 有 效 改善 患 


BH 


疗 使 


且 磁 场 聚焦 和 控制 单个 笔 形 射 束 ， 从 而 允许 在 各 上 


j 笔 形 束 进行 束 流 输送 ， 其 中 扫 


一 个 优势 是 离子 治 
与 光子 相 比 ， 离 子 具有 
第 三 个 优势 是 离子 治疗 中 精确 
描 磁 体 在 与 束 流 方向 正 交 的 方向 上 引导 离子 射 束 ， 


了 利用 布拉格 效应 ， 能 够 将 束 流 能 量 最 大 


更 强 的 生物 效应 外， 在 有 限 


的 剂量 输送 。 离 子 治 


的 能 量 深度 处 逐 层 快速 照 冉 台 区 。 
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子 治 疗 相 比 ， 离 子 治疗 可 以 对 更 复杂 的 肿瘤 结构 产生 精确 的 剂量 输送 。 


当前 离子 束 放射 治疗 计划 通常 使 用 X 射线 计算 机 断层 扫 
将 CT 的 Hounsfield 单位 转换 成 相对 
F X 射线 与 物质 相互 作 
异 ， 使 得 这 种 转化 存在 3.5%-5% 的 误差 。 离 子 CT 是 一 和 有 


肿瘤 上 。 但 由 


阻止 本 领 RSP)， 从 而 使 离子 束 的 Bragg 


直接 应 


通过 


更 加 准确 。 离 子 CT fi 
RAC aa). TR 
制作 治疗 计划 。 由 于 带电 离子 与 物质 相互 作 


确 分 布 四 。 


需要 对 单个 离子 穿 过 患者 的 路 径 进行 估计 ， 
这 通常 需要 两 组 位 置 灵 敏 的 检测 


HA 
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描 CCT) 来 制定 S99， 通过 


峰 落 在 目标 


] 和 离子 束 与 物质 相互 作用 时 发 生 的 物理 过 程 存在 差 
' 绕 过 此 转换 步骤 的 成 像 方法 ， 


它 


离子 作为 成 像 粒 子 并 直接 计算 患者 体内 的 RSP 用 于 剂量 计算 ， 使 得 离子 放疗 
上 共有 足够 初始 能 量 的 高 能 离子 束 穿 过 患者 并 最 终 进 入 测量 离子 残 


测 器 探测 到 的 每 个 离 


che. 


了 的 能 量 损 失 用 于 重建 患 


后 ， 


目前 计算 机 模拟 习 
实现 模拟 。 这 些 模拟 大 多 基于 理论 下 
j 蒙 特 卡 罗 方 法 对 CT 成 像 过 程 的 模拟 成 为 这 一 领域 的 最 佳 选择 。 


本 研究 使 


响 。 


现实 两 种 离子 CT 系统 探测 器 ， 利 


2 ”材料 和 方法 


2.1 


本 研究 使 月 
CT 系统 。 探 测 器 距离 如 
的 距离 为 480mm。 没 有 设置 额外 的 剩余 能 量 
计算 探测 器 D2 和 D3 上 游 表 面 的 离子 能 量 差 来 确 


离子 CT 系统 


即 最 可 能 的 路 径 ， 


的 阻止 本 领 以 用 于 


的 物理 特性 ， 离 子 会 在 物质 中 发 生 散 射 ， 这 


以 实现 能 量 损失 位 置 的 更 准 


器 系统 〈 跟 踪 器 平面 ) ， 分 别 位 于 模 体 的 前 


四 个 位 置 灵敏 探测 器 用 于 测量 最 可 能 路 径 估计 的 离子 入 
F 台 对 CT «jl E et 
究 ， 无 法 真实 的 反映 CT 系统 的 成 像 过 程 。 为 此 ， 利 


| 途径 go 


蒙特 卡 罗 程 序 Geant4/Gate[11][12] 搭 建 了 离子 CT 模拟 平台 ， 设 计 了 理 
最 大 似 然 法 和 ASD-POCS 算法 重建 了 图 人像， 研究 了 理 
想 与 现实 两 种 设置 、 多 种 能 量 下 以 及 不 同 离子 种 类 下 对 模 体 内 四 种 材料 重建 RSP 误差 的 影 


口 和 出 口 位 置 。 
如 利用 Matlab、C++ 等 软件 


想 与 


H Geant4 /Gate 8.2 来 构建 具有 四 个 位 置 灵 敏 探测 器 (D1,D2,D3,D4) 的 离子 


图 1 所 示 ， 两 侧 控 测 器 之 间 的 距离 为 Somm， 内 侧 两 个 探测 器 之 间 


探测 器 ， 直 接 在 位 置 探测 器 上 测量 能 量 ， 通 过 


定 模 体 中 的 


Ak =E. 


能 量 损失 o 


50mm 480mm 50mm 


图 1 蒙特 卡 洛 中 使 用 的 离子 CT 系统 结构 示意 图 

水 等 效 路 径 长 度 (Water Equivalent Path Length， 简 写 为 WEPL) 是 在 离子 成 像 重建 
中 ， 离 子 能 量 损 失 被 转换 为 沿 着 离子 路 径 通 过 模 体 整个 路 径 上 RSP 的 积分 。 探 测 器 对 穿 过 
模 体 之 前 的 初始 能 量 〈Ein) 和 穿 过 模 体 之 后 的 能 量 CE. 进行 评估 。 然 后 根据 公式 DU, 
通过 在 水 中 的 阻止 本 领 上 积分 ， 将 该 能 量 损失 转换 为 WEPL。 


Eout 


wert = | i (D 


gu, (dE/dX)water 

SRIM (Stopping and Range of Ions in Matter) 是 一 个 基于 蒙特 卡 洛 方法 的 计算 程序 ， 
专门 用 于 模拟 计算 离子 在 物质 中 的 阻止 本 领 和 射程 。 通 过 使 用 SRIMII4 中 的 Stopping/Range 
Table 功能 ， 对 0-350 MeV (ER MeV/u) 能 量 范 围 的 质子 和 碳 离子 ， 部 体 材 料 设置 为 密度 
lg/em? 的 水 ， 获 得 积分 所 需 的 阻止 本 领 表 。 

为 了 更 接近 现实 情况 ， 本 研究 使 用 Geant4/Gate 内 置 的 治疗 计划 系统 (TPS) 05 并 参考 
了 在 德国 海德 您 离子 束 治疗 中 心 使 用 的 真实 笔 形 束 扫描 参数 作为 源 09。 详 细 参 数 如 表 1 所 
不 。 


表 1 笔 形 束 设置 


参数 名 称 数值 

光斑 尺寸 (x) (标准 偏差 〉(mm) 3.0 
光斑 尺寸 (y) 《标准 偏差 ) (mm) 3.0 
光斑 发 散 度 0 (mrad) 2.8 
光斑 发 散 度 @ (mrad) 2.8 
光斑 发 射 度 0 (mm*mrad) 3.0 
光斑 发 射 度 o (mm*mrad) 3.0 


22 HM 

本 研究 创建 使 用 了 一 个 圆柱 型 水 模 体 放置 在 原点 ， 直 径 14cm, 12.5cm, HABA 4 
个 直径 lem, f 2.5em 的 圆柱 形 插件 ， 分 别 由 硫 、 磷 、 碳 和 钙 组 成 ， 密 度 分 别 为 2.068 
g/cm?, 1.821 g/cm?, 2.253 g/cm? 和 1.540 gcm3， 如 图 2 所 示 。 


图 2 模 体 示意 图 
通过 使 用 蒙特 卡 洛 计算 程序 SRIM， 利 用 公式 2073， 根 据 不 同情 况 下 离子 束 在 对 应 材 
料 与 水 中 的 阻止 本 领 的 比值 ， 计 算出 不 同 能 量 ， 不 同 离子 种 类 下 4 个 插件 材料 的 相对 阻止 
本 领 (RSP) 值 作为 参考 值 如 表 2 所 示 。 

RSP = -eva (2) 


— YE water 


其 中 ，B 表示 离子 能 量 ，x 表示 离子 在 介质 中 的 路 径 长 度 ，Y 表示 离子 类 型 ，M 表 
示 介 质 材料 种 类 ， 一 (@) 表示 带电 离子 在 介质 中 单位 路 径 长 度 损失 的 能 量 ， 一 CEE)vsw 和 


= 


-Ey r water 分 别 表示 带电 离子 在 材料 和 在 水 中 的 阻止 本 领 。 


表 2 不 同 条 件 下 模 体 中 插件 的 RSP 参考 值 


离子 能 量 (MeV BK MeV/u) 


模 体 材料 离子 种 类 
230 270 310 330 
质子 1.646 1.648 1.650 1.651 
硫 u 
RAT 1.649 
质子 1.418 1.420 1.421 1.422 
BE - 
RAF 1.420 
质子 1.996 1.995 1.994 1.993 
碳 = 
RAF - - - 1.990 
质子 1.214 1.217 1.219 1.220 
E 
碳 离子 - - : 1219 
23 ”数据 采集 


用 离子 束 对 模 体 进行 扫描 ， 模 体 绕 中 心 轴 每 旋转 5° 采 和 集 一 组 数据 ， 每 一 组 数据 使 用 的 
离子 数量 为 993600， 共 获得 72 组 不 同 角 度 投 影 数据 。 
24 ”图 像 重 建 

本 研究 使 用 最 大 似 然 图 像 重 建 方法 03] 重 建 每 个 角度 的 图 像 ， 根 据 其 最 可 能 路 径 估计 将 
每 个 离子 的 WEPL 分 配给 其 通过 的 通道 Kk。 每 个 通道 的 水 等 效 厚度 〈WET) 为 : 


12 
k, 
XN i$ WEPLy, 


WES E S (3) 


每 个 离子 n( 通 过 通道 的 离子 总 数 为 N ) 对 通道 k 的 贡献 具有 加 权 作 用 ， 有 具体 取决 于 
它 在 通过 总 长 度 Ln 时 在 某 个 通道 中 花费 的 长 度 1。Ln 为 两 个 最 内 层 探测 器 平面 之 间 的 距 
离 。 由 此 ， 得 到 每 个 角度 的 图 像 数据 。 
使 用 得 到 的 每 个 角度 的 图 像 数 据 作为 输入 ， 参 考 Kaser 等 人 的 工作 设置 ， 利 用 ASD- 
POCS 算法 HI， 和 迭代 次 数 设 为 15 次 ， 以 1xlxlmm 重建 CT 图 像 。 

2.5 评估 
使 用 模 体 的 RSP 重建 图 像 分 析 不 同 材料 的 RSP 重建 值 与 参考 值 的 误差 09。 选 取 每 个 
插件 中 间 切 片 6x6 C6 mmx6 mm) 像素 区 域内 所 有 像素 RSP 平均 值 为 RSP 重建 值 ， 以 相对 
误差 表示 。 公 式 如 下 : 


— 
ni 
in 


RSP aaa 7 RSP 参考 值 | 


| 
相对 误差 = x 100% 


3 结果 


3.1 ”探测 器 对 RSP 重建 误差 的 影响 

为 了 探究 研究 了 探测 器 对 RSP 重建 误差 的 影响 ， 设 置 了 两 种 探测 器 配置 。 第 一 种 为 理 
想 设 置 ， 即 为 了 尽量 减少 探测 器 材料 对 重建 结果 的 影响 ， 将 跟踪 探测 器 的 材料 设 为 空气 ， 
从 而 不 对 离子 的 能 量 和 位 置 方 向 造成 影响 。 第 二 种 为 现实 设置 ， 即 为 了 研究 在 更 现实 条 件 
下 的 重建 性 能 ， 将 每 个 跟踪 探测 器 材料 为 硅 ， 厚 度 设 置 为 双 面 硅 微 条 探测 器 典型 尺寸 的 
300 hm。 两 种 设置 的 唯一 区 别 仅 为 探测 器 材料 不 同 。 
使 用 能 量 为 330MeV 的 质子 分 别 在 两 种 设置 的 系统 下 对 模 体 进行 照射 重建 ，RSP E 
建 图 如 图 3 所 示 ， 四 种 插件 材料 的 RSP 相对 误差 如 图 4 所 示 。 理 想 设置 和 现实 设置 下 硫 、 
磷 、 碳 、 钙 四 种 材料 的 相对 误差 分 别 是 0.181 %、0.351%、0.250%、0.245% 和 和 0.787%. 
1.547%. 0.65296. 0.65596. 
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b， 现 实 设 置 RSP 重建 区 
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a. 理想 设置 RSP 重建 


图 3 两 种 设置 RSP 重建 图 
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图 4 两 种 设置 的 RSP 相对 误差 
由 图 4 可知， 两 种 设置 下 各 材料 的 RSP 相对 误差 均 小 于 1.547%， 理 想 设 置 下 的 RSP 
相对 误差 远 小 于 现实 设置 的 误差 ， 在 钙 材料 中 的 RSP 相对 误差 为 0.245%， 相 同 条 件 下 现 
实 设置 的 误差 达到 0.655%， 接 近 理想 设置 误差 的 三 倍 多 。 为 了 排除 探测 器 对 离子 束 流 的 干 
扰 并 影响 成 像 结果 ， 后 续 的 研究 红 使 用 理想 设置 进行 。 
3.2 ”不 同 能 量 离子 对 RSP 重建 误差 的 影响 

分 别 使 用 能 量 为 230MeV. 270MeV. 310MeV 和 330MeV 的 质子 对 模 体 进行 照射 重 
Æ, RSP 重建 图 如 图 $ 所 示 ， 四 种 插件 材料 的 RSP 相对 误差 如 图 6 所 示 。 
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图 6 不 同 能 量 RSP 相对 误差 

由 图 6 可 知 ， 各 材料 重建 RSP 的 相对 误差 随 着 质子 能 量 升 高 而 减 小 ， 在 质子 能 量 大 
230MeV 时 最 大 ， 材 料 硫 、 磷 、 碳 、 钙 中 的 重建 RSP 相对 误差 分 别 为 2.855%、2.468%、 
1.653%、2.553%， 在 330MeV 时 的 RSP 相对 误差 达到 最 小 ， 在 硫 、 磷 、 碳 、 钙 四 种 材料 中 
分 别 为 0.181%、0.351%、0.250% 和 0.245%. 
3.3 ”不同 离子 种 类 对 RSP 重建 误差 的 影响 

分 别 使 用 目前 离子 放疗 中 常用 的 两 种 离子 为 质子 与 碳 离 子 ， 以 能 量 为 330MeV (ak 
MeV/u) 对 模 体 进行 照射 重建 ，RSP 重建 图 如 图 7 所 示 ， 四 种 插件 材料 的 RSP 相对 误差 如 
图 8 所 示 。 碳 离子 和 质子 在 磷 、 磷 、 碳 和 钙 材 料 中 的 RSP 相对 误差 分 别 为 0.060%、 
0.281%、0.150%、0.082% 和 0.181%, 0.351%, 0.250%, 0.245%. 
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7 不 同 离子 种 类 RSP 重建 图 
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图 8 不 同 离子 种 类 RSP 相对 误差 
由 图 8 可 知 ， 碳 离子 束 在 330 MeV/u 能 量 时 的 RSP 相对 误差 都 在 0.281% 之 内 ， 远 小 
于 相同 条 件 下 质子 的 RSP 相对 误差 ， 在 材料 钙 中 获得 的 重建 RSP 相对 误差 为 0.082%， 在 
硫 中 的 RSP 相对 误差 为 0.060%， 相 同 条 件 下 质子 束 在 钙 和 硫 中 重建 得 到 的 RSP 相对 误差 
分 别 为 0.245% 和 0.18196. 


4 讨论 

本 文 基于 蒙特 卡 罗 模 拟 研究 了 离子 束 CT 在 不 同 条 件 〈 理 想 与 现实 离子 CT 系统 ， 不 
同 能 量 ， 不 同 离子 种 类 ) 对 图 像 重 建 及 RSP 误差 的 影响 。 在 离子 CT 系统 探测 器 对 重建 
RSP 误差 影响 的 研究 中 ， 现 实 条 件 下 硫 的 RSP 相对 误差 为 0.787%， 理 想 条 件 下 为 
0.181%， 有 着 四 倍 多 的 误差 ， 两 者 设置 在 其 他 材料 的 相对 误差 差距 也 很 明显 ， 说 明 在 离子 
CT 的 现实 建造 中 ， 如 何 减 小 离子 CT 系统 探测 器 对 于 离子 束 的 库仑 散射 影响 从 而 得 到 更 加 
准确 的 RSP 图 像 将 成 为 关键 。 本 研究 也 存在 着 一 些 局 限 性 ， 模 拟 与 现实 存在 差异 ， 后 续 可 
将 探测 器 响应 和 扫描 过 程 中 的 辐射 环境 纳入 考虑 。 


5 结论 


本 研究 通过 使 用 蒙特 卡 罗 程 序 Geant4/Gate 搭建 了 离子 CT 模拟 平台 ， 设 计 了 理想 与 现 
实 两 种 离子 CT 系统 探测 器 ， 利 用 最 大 似 然 法 和 ASD-POCS 算法 实现 了 图 像 重 建 ， 研 究 了 
理想 与 现实 两 种 设置 、 多 种 能 量 下 以 及 不 同 离子 种 类 下 对 模 体 重建 RSP 误差 的 影响 。 得 出 
以 下 结论 : (1) 330MeV 质子 在 理想 与 现实 设置 两 种 情况 下 的 RSP 相对 误差 都 小 于 
1.547%， 理 想 设 置 下 的 RSP 重建 误差 远 小 于 现实 设置 ， 现 实 设置 下 各 材料 的 RSP 重建 误 
差 接近 理想 设置 下 的 三 倍 ; (2) 质子 的 RSP 相对 误差 随 着 能 量 的 增加 而 减 小 ， 在 
230MeV 时 的 RSP 相对 误差 最 大 ， 在 硫 、 磷 、 碳 和 钙 四 种 材料 中 分 别 为 2.855%、2.468%、 


1.653%, 2.553%. TE 330MeV 时 的 RSP 相对 误差 达到 最 小 ， 在 硫 、 磷 、 碳 和 钙 四 种 材料 中 
分 别 为 0.18196. 0.35196. 0.250%Fll 0.245%; (3) 在 能 量 330MeVA 下 ， 碳 离子 在 硫 、 
磷 、 碳 和 钙 四 种 材料 中 的 RSP 相对 误差 分 别 为 0.060%、0.281%、0.150%、0.082%， 误 差 
均 在 0.281% 以 内 ，RSP 相对 误差 小 于 能 量 330MeV 的 质子 。 因 此 ， 碳 离子 CT 应 比 质子 
CT 更 能 为 离子 放疗 治疗 计划 设计 提供 准确 的 RSP 数据 。 
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Study on ion CT image reconstruction and 
relative stopping power error based on Monte 


Carlo simulations 
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Abstract: Ion beam radiotherapy is increasingly being used for tumor treatment. Currently, the 
use of traditional X-ray computed tomography (CT) for treatment planning of ion-beam 
radiotherapy has significant relative stopping power (RSP) errors. A more ideal approach is to 
directly use ion beams to generate patient images for treatment planning so as to avoid RSP 
conversion and reduce RSP errors. In this study, the Monte Carlo program Geent4/Gate was 
adopted to establish an ion-CT simulation platform, designed two sets of ion CT systems, ideal 
and real, and reconstructed images using the maximum likelihood method and ASD-POCS 
algorithm. The effects of the ideal and real settings, multiple energies, and different ion types on 
the RSP error of phantom reconstruction were investigated. The results show that the relative error 
of RSP for 330MeV protons in both the ideal and real settings was less than 1.547%, and the RSP 
reconstruction error in the ideal settings was much smaller than that in the real settings. The RSP 
reconstruction error of each material under realistic settings is close to three times that under ideal 
settings. The relative errors of RSP for protons decrease with the increase of incident energy. The 
relative errors of RSP were the biggest at 230MeV, and were2.855%, 2.468%, 1.65396. and 
2.553% in sulfur, phosphorus, carbon, and calcium materials. The RSP relative error reached its 
minimum at 330MeV, with 0.181%, 0.351%, 0.250%, and 0.245% in sulfur, phosphorus, carbon, 
and calcium materials. At energy of 330MeV/u for carbon ions , the RSP relative errors in sulfur, 
phosphorus, carbon, and calcium materials were 0.060%, 0.281%, 0.150%, and 0.082%, 
respectively, all within 0.281%. And the RSP relative errors were much smaller than those of 
protons under the same conditions. Thus, compared with proton CT, carbon ion-beam CT seems 
more possible to provide accurate RSP data for treatment planning for ion-beam radiotherapy. 
Keywords: Ion CT; Monte Carlo simulation; Proton; Carbon ion; Relative stopping power 
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